a. I 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2887 


ON THE STABILITY OF THE LAMINAR MIXING REGION 
BETWEEN TWO PARALLEL STREAMS IN A GAS 

By C. C. Lin 

Massachusetts Institute of Technology 



Washington 
January 1953 





O-' 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


I i-i/ii LlsnAiu ru\i a, NM 



00bS7^S 


TECHNICAL NOTE 2887 

ON THE STABILITY OF THE LAMINAR MIXING REGION 
BETWEEN TWO PARALLEL STREAMS IN A GAS 
By C. C. Lin 

SUMMARY 


The stability of the mixing of two streams was studied both for 
the interest in the problem and for clarifying certain points in the 
basic stability theory. It is shown that, when the relative speed of 
the two parallel streams exceeds the sum of their velocities of sound, 
subsonic oscillations c ann ot occur and the mixing region may be expected 
to be stable with respect to small disturbances. It is further shown 
that, when viscosity and heat conductivity are neglected, if the flow 
ran execute a small neutral subsonic oscillation of finite, wave length, 
it can also execute self -excited oscillations of longer wave lengths 
nnd damped oscillations of shorter wave lengths. 

Rigorous developments of the mathematical theory of asymptotic 
solutions confirm previous methods of solution of the stability equa- 
tions in a compressible fluid. This theory also shows that, at high 
Reynolds numbers, the damped osci ll ations in a strictly parallel main 
flow have a structure similar to that of the vorticity field in fully 
developed turbulent flow. 

Sample calculations are also included exhibiting various quantita- 
tive properties of these small oscillations. 


INTRODUCTION 


The mixing of two parallel streams of gas occurs in a number of 
cases. An interesting example is furnished by the slip stream in a 
three-shock configuration. It has long been suggested that such la min a r 
rrH vtng zones could, at sufficiently high speeds, be stable with respect 
to small disturbances although they are known to be very unstable at low 
speeds. The purpose of the present investigation is to find out some 
of these stability characteristics. Apart from the development of the 
general theory, there are included the calculations of the neutral and 
unstable oscillations, the extent of the amplification, and other related 
properties . 
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The basic equations for the study of small disturbances in the 
laminar boundary layer of gases have been given in reference 1. However, 
the development of the theory there has specific reference to the case 
of a layer near a solid bo undar y. As pointed out in references 2 and 3> 
the stability theory for a mixing zone in an incompressible fluid dif- 
fers from that for a layer near a solid boundary in that solutions of 
the exponential type are unimportant. This leads to the conclusion that 
the effects of viscosity and heat conduction are negligible, except at 
very low Reynolds numbers, in determin in g the characteristics of the 
oscillations. To’ confirm this point, a rigorous mathematical theory of 
asymptotic solutions was developed for the compressible case similar to 
that indicated in reference k for the incompressible case. In view of 
the mathematical interest involved, it was decided that this basic part 
of the present investigation would be published separately in mathematical 
journals (see references 5 and 6), and only the main results and their 
physical significance will be presented here. 

As in the case of an incompressible fluid, the "inviscid" case is 
expected to be characteristic of the behavior of the disturbances at 
moderately large Reynolds numbers. Most of the studies are, therefore, 
made in the inviscid case. However, the interpretation of the inviscid 
case must be subjected to the same care as in the incompressible casej 
that is, in the case of damped disturbances, the differential equation 
of the inviscid flow may not be regarded as valid throughout the real 
axis. There is a finite viscous region even in the limit of vanis hin g 
viscosity. The complex conjugate of the amplified disturbance is cer- 
tainly a solution of the inviscid equation, but it is not a limiting 
solution of the complete viscous equation. This behavior of the inviscid 
solution reminds one of the vorticity structure of fully developed tur- 
bulent flow as found by Batchelor and Townsend (reference 7) . 

The nonexistence of subsonic- disturbances is visually associated 
with the stability of the parallel flow. There seems to be some basis 
for doing this, although the role of supersonic disturbances has never 
been fully clarified. It Is easy to see that, for certain combinations 
of the properties of the two streams, it is impossible to have a sub- 
sonic disturbance relative to both. Under such conditions, one may 
expect stability. These conditions for stability are developed herein 
and are expected to hold, irrespective of the viscous effects. 

Applying the theory of stability in the inviscid case, one can 
further narrow down the possible range of instability. This will depend 
upon the velocity and temperature distributions in the shear zone. In 
the present work, calculations are made for gases with Prandtl number 
equal to unity. Although the condition of equal total enthalpy in the 
two streams is also used, it is shown that this restriction can be 
immediately removed by considering a moving observer. It is found 
that the condition of stability thus found does not differ very much 
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from that found above from general considerations. Thus, it may he 
surmised that the exact distributions of temperature and velocity have 
only a secondary influence on the stability characteristics in the 
mixing zone. Thus, the approximations used in the present calculations 
of the basic velocity and temperature distributions cannot influence 
the final results to any appreciable extent. 

The neutral disturbances are of two kinds: (l) A steady deviation 

and (2) an osci ll ation of finite wave length. Thus, there are two 
branches of the neutral curve at infinite Reynolds numbers . They may 
be expected to join together at low Reynolds numbers enclosing a region 
of instability. Calculations of neutral and amplified oscillations are 
carried out in a number of cases with one stream at rest. 

This investigation, carried out at the Ifessachusetts Institute of 
Technology, was sponsored by and conducted with the financial assistance 
of the Rational Advisory Committee for Aeronautics. The author is 
indebted to Mr. D. W. Dunn for his valuable suggestions and help in the 
preparation of the final version of the report and to Miss Diana Mason 
and Mr. W. V. Caldwell for their help in making the numerical calculations. 


STEADY FLOW IN THE LAMINAR LAYER BETWEEN 
TWO PARALLEL STREAMS 


The basic steady flow under discussion is a boundary- layer flow with 
no body forces and no pressure gradient. The basic equations are (see 
list of symbols in the appendix) : 
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the pressure "being a constant throughout the field. In the case of a 
homogeneous incompressible fluid, with the streams at the same tempera- 
ture, the temperature may he taken as constant throughout the field, 
and the first two equations can be solved for the velocity distribu- 
tions u*(x*,y*) and v*(x*,y*) . In the case of a compressible gas, 
the integration has to be carried out for in dividual cases. However, 
if the Prandtl number Cp|i*/k* is equal to unity, it is known that 

there is a quadratic relation of the type 


CpT* + i u * 2 = A + Bu* 


(k) 


between the temperature and the velocity, and one is again essentially 
dealing with two distributions u*(x*,y*) and v*(x*,y*). Indeed, the 
constants A and B are given in terms of the conditions in the _ 
parallel streams as follows: 


B - |( u 2 + u l) + =p( T 2 - T l)/( u 2 - %) 

A = - | U X U 2 + - T 2 0 1 )/(U 2 - Uj 


(5) 


If the total enthalpy in the two streams 1 and 2 is the same, that is, 
if 


OpTl + | V - CpT 2 + i U 2 2 


( 6 ) 


then B = 0, and the total enthalpy is constant throughout the whole field: 

OpT* + | u * 2 = C 


(7) 
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This is a particularly simple case, to which, however, all other cases 
can he easily reduced. This is done hy rewriting equation (k) into 
the form 


CpT* + i(u* - B) 2 = A + B 2 /2 
2 


( 8 ) 


One need only consider an observer moving with the speed B and con- 
sider the relative velocity u* - B. Thus, the solution for general 
specified values of Up,T^, and U 2 ,T 2 can he derived from the iso- 

energetic solution with boundary conditions Up - B,Tp and U 2 - B,T 2 
hy simply adding the constant B to the u*- component of the velocity. ^ 

In the following discussions, isoenergetic basic solutions will he 
referred to often; however, it should he kept in mind that hy the con- 
sideration of a moving observer the general case may he obtained." This 
transformation is not restricted to the steady flow hut applies to the 
consideration of the disturbances as well. Thus, if all the cases of 
constant total enthalpy are calculated, all the other cases are also 
known. 

So far, the viscosity coefficient may depend on the absolute tem- 
perature In any manner. If there is direct proportionality of these 
quantities, the solution in the compressible case can he expressed in 
terms of that in the incompressible case. These relations are well- 
know and, in the following discussion, only the results relevant to 
this case will he given. 


Incompressible Case 
For the incompressible case: 


u* = U-jf '( ti) 

v* = f^f' - f) 


(9) 


lit is to be noticed that in certain cases with T 2 /T;i_ < 1 the 
boundary value U2 - B for the corresponding isoenergetic problem may 
be negative while Up - B is still positive. Thus the isoenergetic 
problem may not be physically significant. However, for the purposes 
of theoretical analysis, this point is not important. 
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where 



and’ f( tj) satisfies the differential equation 


f f ' ' + 2f ’ 1 ’ = 0 


( 11 ) 


with the conditions 


f * (Tj) -> 1 as T] > oo 


^ as tj 
U 1 


( 12 ) 


A third condition is arbitrary tip to a translation along the Tj-axis . 

For a typical scale, the momentum-boundary- layer thickness 0* may 
he introduced, which is defined hy 


/ OO 

p*(u x - u*)(u* - U 2 ) dy* 


(13) 


Then 


0 = 




f’(Tj) - ^ 


u 2 


•Uo 


dTj 


(14) 


for an incompressible fluid. It will he seen later that the same for- 
mulas apply for a compressible fluid in the isoenergetic case. 

Calculations made hy Gortler (reference 8) for the turbulent mixing 
region can he easily adopted for the purpose at hand. There is only a 
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slight difference in the method of representation. To convert his 
function F'(|) to the notation of this report, the following rela- 
tions should he used 


f’(tl) = -SL-. p*(|) 
1 + X 

t] = 2jmi 

where 



Ui +u 2 

The converted results are given in table I and figure 1. 


(15a) 


(15b) 


Compressible Case 

For the compressible cane with constant total enthalpy, p* cc T*, 

where 


u* = Ujf’U) 



(16) 


(17) 


(18) 


where 


t\ = v^/u-l 
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The momentum thickness 9 * can "be defined by the same equation (equa- 
tion ( 13 )) and can be shown to have the same value. (equation (l4) ) in 
dimensionless form. Velocity and temperature distributions in the case 
U2 = 0 have been calculated for several Mach numbers) these are tabu- 
lated in table II and plotted in figure 2. 

A GENERAL RESTRICTION FOR EXISTENCE OF 
NEUTRAL SUBSONIC DISTURBANCES 


Before going into the general theory of stability, a preliminary 
discussion of the stability of the mixing zone will be given. It will 
be found below that the results thus obtained give quite an adequate 
description of the general dependence of the stability of the mixing 
zone on the Lfe.ch numbers of the streams. It will be shown that, if 
the average Mach number "M of the relative motion (defined by equa- 
tions (25) and (26)) exceeds the value 2, subsonic disturbances in the 
sense of reference 1 cannot exist, and the mixing zone may be expected 
to be stable. 

Consider two parallel streams at speeds U^ and U 2 , tempera- 
tures T^ emd T 2 , and Mach numbers Mj_ and M 2 . For definiteness, 
take - U 2 > 0. Consider an observer moving with the speed U 2 . 

Then the streams appear to have speeds Up - U 2 and 0, while their 
temperatures are obviously not changed. If 


c' = c - U 2 


(19) 


denotes the speed (relative to the moving observer) of a wavy motion 
propagating in the direction of the stream (c 1 positive or negative), 
then the conditions for subsonic disturbances are 


|( u l - W 2 ) - o' | < aj. 


( 20 ) 


| c ' j < a 2 (21) 

where and a 2 are the velocities of sound in the streams. It is 

clear that subsonic disturbances can always exist if U-j_ - Ug is less 
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than a-]_. For - Ug > a-j_, condition (20) can he satisfied only 
with c' > 0. Thus 

0 < c' < a 2 , that is, 0 < c - U 2 < a 2 (22) 

Now, note that condition (20) is actually, 

|Ul - c| < a i 

and hence 

_ - &1 < U x - c < a 1 (23) 

Adding relations (22) and (23), it is found that 


U 1 “ U 2 < a l + a 2 


(210 


is a necessary condition for the existence of subsonic disturbances. 
By introducing the average velocity of sound 


(25) 



and the average Mach number of relative motion 


M = 



(26) 


defined with respect to this average sound velocity, it may be concluded 
that subsonic disturbances cannot exist (and the mixing zone may be 
expected to be stable) if the average Mach number M of the relative 
motion exceeds the value 2, as stated at the beginning of this section. 
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In terms of the Mach numbers of the two streams, condition (24) 
states that if 


Ml > 1 + \|t 2 /Ti(i + Mg) 


( 27 ) 


Subsonic disturbances cannot exist, and the motion is stable. (See 

fig. 3 .) 

In the case of constant total enthalpy, the temperature ratio 
Tg^T^ can expressed as 


Tg _ 1 + (7 - 1 )Mx 2 /2 

Tl 1 + (7 - l)Mg 2 /2 


Then condition (27) becomes 


Ml >1 + 


1 + (7 - 1)M x 2 /2 
1 + (7 - l)M2 2 /2 


11/2 


( 1 + M 2 ) 


The critical condition is 


( 28 ) 


(29) 


M-l = 1 + 


1 + (7 1)^/2 

1 + (7 - l)Mg 2 /2 


1/2 


(1 



(30) 


Removing the square root. 


1 + (7 - 1)M-L 2 j 

,2 1 

[Mi - l) 

1 + (7 - l)M 2 2 j 

1 2 * 

[“2 +1 ) 


This form suggests an obvious solution, Mp = -Mg . Another solution 
can then be easily obtained as 


Ml = 


(3 - 7)Mg + 4 
(3 - 7) - 2(7 - l)Mg 


(31) 
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By substituting in equation ( 30 ), Mp = -Mg is found to "be an extra- 
neous solution. The only solution of equation ( 30 ) is then given by 
equation ( 31 ) . 

Thus the flow is stable if 


(3 - 7 )M 2 + 4 

Mp > 

(3 - 7 ) - 2(7 - l)Mg 


(32) 


The curve for equation ( 31 ) is shown in figure 4, where the region to 
the left of curve A is the region of possible instability. Note the 
symmetry with respect to the line Mp = -Mg, corresponding to a change 

of the positive direction of the x-axis in the physical problem. The 
line Mp = Mg is drawn to take care of the condition Up - Ug > 0, 

since Up = Ug when Mp = Mg. The asymptotes of the curve for equa- 
tion ( 31 ) are 


Mg = (3 - 7 )/ 2(7 - 1) 


and 


Mp = -(3 7)/2(7 - I) 


With 7 = 1.1)., these values are ±2. The intercepts are 


Mp = V(3 - 7) 


and 


Mg = -4/(3 - 7) 


(33) 


(34) 


(35) 

(36) 


With 7 = 1.4, these values are ±2.5. (Cf. reference 9 .) 
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GENERAL STUDY OF SMALL DISTURBANCES IN A NEARLY PARALLEL 
FLOW F IELD IN A COMPRESSIBLE FLUID 


The general theory of a small disturbance in a field of nearly par- 
allel flow of a gas has been developed in reference 1. The rigorous 
mathematical proof and improvement of the theory are given in detail in 
references 5 and 6. In this report, merely the main conclusions and 
their physical interpretations are outlined without going into the 
details. Applications of the theory to the specific case at hand will 
be discussed in some detail. 

Consider a nearly parallel stream with dimensionless velocity and 
temperature distributions w(y) and T(y) . The neglect of the depend- 
ence of these quantities on x and the omission of the y-component of 
the basic flow can be justified by detailed investigations. The line- 
arized differential equations for small disturbances then possess solu- 
tions of the type 


u' 

v' 

P' 

P’ 


Re|f(y)e :kt '( x - Ct il 

Re |^(y)e la ( x - ct ]] 
Re [r(y)e^ a ( x " c ^^] 

Ref(y)e la '( x - Ct i] 


T 1 


Kef(y)e ia ( x - Ct 3 


(37) 


where u', v', p', p', and T l are the perturbations of the two com- 
ponents of velocity, the density, the pressure, and the temperature, add- 
in a suitably defined dimensionless form. The constants a and c are, 
respectively, the real wave number and the complex wave speed. 

The differential equations for the amplitude functions f(y), 0(y) , 
r(y), jt(y), and ©(y) are rather complicated. However, if the effects 
of viscosity and heat conduction are neglected, they become a relatively 
simple system, which can be reduced to the following single differential 
equation for 0(y) : 
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(w - c) 0 ' - w '0 


d 

<3y|T - Mj_ 2 (v - c)^J 


a 2 (w - c )0 


T 


= 0 


(38) 


All the other variables can be expressed in terms of 0 and 0’ as 
follows : 


f = -i[m 2 (w - c)w'0 - T0' I T - M 2 (w - 0)3 
r = ij^ 0 ' + if) + p' 0 j^(w - c) 
it = i 7 M 2 p|^i(w - c)f + w' 0 j 


> 


e = t(Vp - r /p) 


(39) 


The boundary conditions are that the disturbance should be bounded 
as y 7 ±<n. 

The inviscid system would have given a well-defined characteristic- 
value problem if it were not for the fact that differential equation (38) 
has a singularity at the point y = y c where w(y) = c. This singu- 
larity disappears only when d(pw') /dj - 0 at the same point. Other- 
wise, a solution of equation (38) has a logarithmic singularity at 
y = y c and the characteristic-value problem associated with this equa- 
tion becomes indeterminate until the proper branch of the solution is 
determined . 

The determination of the proper branch of the solution and its 
associated physical interpretation is one of the most delicate points 
in the theory of hydrodynamic stability. The mathematical analysis of 
the solutions of the complete viscous equations and their limiting solu- 
tions will be made first before discussing their physical interpretation. 

The complete system of viscous equations can be shown to be equiva- 
lent to a system of six linear equations of the first order in six 
unknowns. Thus, there are six independent solutions. These solutions 
have been formally obtained as asymptotic series in reference 1 , and 
their rigorous mathematical investigation has been carried out in 
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references 5 and 6. It is found that two of the six solutions can "be 
expressed in asymptotic series of the form 


f = f(°) + — f^ 1 ) + . . . 
X 2 


0 = 0 (°) + — 0 ( x ) + . . . 

X 2 


( 1 + 0 ) 


where f(^), 0(0), . . . are the inviscid solutions satisfying equa- 

tions (38) and (39) and X 2 = aR, where R is the Reynolds number based 
on the thickness of the mixing region. Thus, the formal limit of equa- 
tions (40) does approach the inviscid solution, hut a complete study of 
these equations also carries the knowledge of the proper branch to be 
used. 


Four other solutions of the complete system of viscous equations 
are of the form 


where 


f = F exp^XQj^ 
0 = <5 exp(XQ i ) 


i = 1, 2, 3, 1+ V 


( 1 + 1 ) 


Qi = -Q 2 






(w - c) dy 


> 


(1+2) 


and F and <J> can be expressed as power series of X, involving only 
a finite number of positive powers. 
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In the establishment of these asymptotic solutions, it is shown that 
the lines Re^Q-^ = 0 , Re^Q^) = and Re^Q-^ - = 0 are of interest. 

The geometry of these lines relative to the point y c (where w = c) and 

the real axis of the y-plane are shown in figure 5. There are asymptotic 
solutions expressed in equations (40) and (4-1) which maintain the same 
analytical expression on the two sides of the dotted lines. However, in 
crossing the solid lines, they generally change their behavior. Thus, 
the following conclusion may be drawn: The proper branch of the multiple- 

valued asymptotic solutions is obtained by taking a path in the complex 
y-plane below the point y = y c (in this case w(y) is monotonically 

increasing along the real axis) . 

This is the branch taken in reference 1. By examining the behavior 
of solutions of the type given in equations ( 4 l), it can be easily shown 
that they diverge for either positive or negative large values of y. 

These solutions should, therefore, be rejected in the present problem. 

The effect of viscosity is then to be obtained through solutions of the 
form of equations (40) . 


GENERAL STABILITY CHARACTERISTICS IN 
THE INVISCID CASE 


Much of the discussion of the inviscid case in reference 1 applies 
to the present case. However, as noted above, in the present problem, 
the possibility can be more readily realized that subsonic disturbances 
may not exist at all, and that the motion may then be expected to be 
completely stable with respect to small disturbances. Another main 
difference lies in the "steady" deviation, that is, solutions of equa- 
tion ( 38 ) with a = 0. In the present discussion, the general line of 
discussion in reference 1 will be followed. However, the difference in 
the boundary conditions often causes a difference in the method of 
analysis. The arguments are, therefore, presented in some detail. Also, 
it will be shown that disturbances having wave-lengths slightly longer 
than that of the neutral subsonic disturbances of finite wave length 
are unstable while those with slightly shorter wave lengths are the 
stable ones. This conclusion also applies to the case of the boundary 
layer, but it was not obtained in reference 1. The analysis also leads 
to an approximate estimation of the dependence of amplification on wave 
length. This will be used for the calcu l ation of the anqplification of 
the disturbances in the section "Self -Excited Oscillations." 

Some general analytical properties of the solutions of equation ( 38 ) 
will be first summarized, particularly for the case of subsonic disturb- 
ances . In this case. 
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T - M]_ 2 (w - c) 2 > 0 (43) 

for both free streams. For large positive values of y, equation (38) 
may he approximated hy 


0' ' - p x 2 0 = 0 with = a 2 1 - M x 2 (l - c) 2 


(44) 


For large negative values of y, it may he approximated hy 


0" - p 2 2 0 = 0 with p 2 2 


a 




(45) 


Since hoth p-^ and p 2 2 are positive, the solution 0(y) is exponen- 
tial in nature for large values of y in the case of subsonic disturb- 
ances, with the exception of the case a = 0. In that case, two inde- 
pendent solutions of equation (38) are 


*= w - c 


(46) 


and 


4 


2 


(w - 



(^ 7 ) 


The first solution is hounded while the second varies linearly with y 
for large values of y. 

For any value of c, the hounded solution (equation (46)) corre- 
sponds to no disturbance at infinity. In fact, for a = 0, the disturb- 
ance v' is identically zero, hy equation (37) • The other components 
of the disturbance are given hy equations (39)* and it can he easily 
verified that they all vanish at large distances. 

It can he shown from the general nature of the temperature and 
velocity distributions that condition (43) holds throughout the mixing 
zone if it holds in hoth free streams. It is then obvious that 
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equation ( 38 ) has a singularity only at w - c = 0. Actual l y, even 
if T - M-j^Cv - c ) 2 = 0 , it only gives rise to an apparent singularity. 

Analysis of the solution in thie neighbor hood of the point y = y c 
(where w = c) gives the following two solutions: 


0 i = (y - y c )si(y - y c ) 


02 = 82 (y - y c ) + K 0 i lQ Se (y - y c ) 


where 



(>tS) 


(49) 


and g 1 and g 2 are power series in (y - y c ) with g-^O) = w c ' £ 0 

nnrl g 2 (o) = T c/ W c ' 0 °* 1116 Proper branch of the logarit hm ic function 

is to be taken in accordance with the method discussed in the preceding 
section. 

For . real values of c, it can be shown that the Reynolds shear 
stress - 


t = -pu’v' 


(50) 


is a constant except for a possible jump at y = y c - fa fact, this 
jump is 

[ T ] =!«(v) 2 K,| 2 /t c 2 

when one passes from y c - 0 to y c + 0. Thus, for a neutral disturb- 
ance with a £ 0 , the condition 

K = 0 (52) 


must be satisfied, since t = 0 for y — ±» and, therefore, does not 
have any jump. 
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Thus, for the existence of a subsonic disturbance, the quantity. 



(53) 


must vanish at some point in the field; furthermore, the corresponding 
value w s of w must he subsonic relative to both streams; that is. 


1 - 


Mi 


< 



(5k) 


The above reasoning can be applied to the case of the boundary layer; 
a somewhat different argument was used in reference 1. A little' calcu- 
lation will show that, in general, condition ( 53 ) will be satisfied, 
although condition (54) may not. If the latter is also satisfied, then 
there actually exists a subsonic disturbance with c = w 8 . To prove 
this, equation ( 38 ) is rewritten in the form 


where 



(55) 


h(y) = [t - M x 2 (w - og - 1 >0 



(56) 


With c = w B , q(y) is regular all along the real axis. This is a 

2 

characteristic-value problem for the parameter k = -a and is associ- 
ated with the variational principle 



dy = 0 


(57) 
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with 



Y dy = Constant 


(58) 


The least value for -a 1 2 is then given hy the minimum value of the 
ratio 



(59) 


for al 1 functions f(y) such that the integrals in this ratio are con- 
vergent. Wow, it is necessary only to show that for certain functions 
the ratio of the integrals is negative. 

In the choice of such a function, it is necessary for only the 
numerator in equation (59) to he convergent. The convergence of the 
integral in the denominator is immaterial, for one can always modify 
the function at sufficiently large values of |y| so that the denomi- 
nator becomes convergent without altering the sign of the numerator. 

To choose a function so that the integral 


r oo 

(hf ,2 + q.x ) dy ( 60 ) 


is negative, it is first noted that 


1 = 0 for f = w - c ( 6 l) 

This follows from the fact that 0 = w - c is a solution of equation (55) 
when a = 0. It can also he directly verified. Obviously, the value of 
I is not changed if f is now changed to |w - cl . A further modifi- 
cation of the function would yield the desired result. In the neighbor- 
hood of y c , | w — c | is small, hut w' is finite. If in this neighbor- 
hood the function f(y) = |w - cj is replaced by a horizontal straight 
line, the integrand in equation ( 60 ) is certainly decreased (since | >0 
nnfl a change in f ' causes a larger change in the integral than that 
in f), and the resultant integral is negative. This completes the proof 
desired. 
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Thus, the sufficiency of conditions ( 53 ) and (54) for the existence 
of subsonic disturbances is also established. Let the corresponding 
value of a he denoted by <Xg . 

Next, it will be shown that for a slightly less than Og the 
disturbances are unstable. To prove this the quantity (dc/dk) s will 
be calculated, anri its imaginary part will be shown to be positive. 

Consider a characteristic function 0(y;k,c) of equation (55) 
corresponding to a given value of k. As c changes, k also changes. 
If equation (55) is differentiated with respect to k, the following 
equation is obtained 


_d_ 

ay 



*£ + t = 

dk T 


0 


(62) 


where 

m i = M + M *£ (63) 

8k 8c dk 


Now, multiply equation ( 62 ) by 0, equation (55) by 0, subtract the 
results, and then integrate with respect to y along a path in the com- 
plex plane which leads from y = -«> on the real axis to y = °° on the 
real axis but passes below the point y c . The following equation is then 
obtained for dc/ dk: 


ic P 8h/d0\ 2 + 8q gj 2 
dk J ^c\dy/ 8c _ 

So far equation (64) is general. Now, specialize to the case of 
the neutral disturbance in question. Most of the integration can then 
be carried out along the real axis, with real resultant values. However, 
the integral 


ay 


-/ 


r 

T 


ay 


(64) 
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must be evaluated along an indented patb with a small detour below y , 

since 8q/8c has a pole at that point. Calculations show that the 
imaginary part of J is 


Im(j) 



( 66 ) 


if this q uan tity does not vanish. Indeed, this quantity can be seen 
to be negative. Thus, equation (64) yields a relation of the form 


— (P + iQ) = R > 0 

dk 


where Q = Im(j) . Thus, 


dc 

dk 


R 

P 2 + q2 


(P - iQ) 


has a positive imaginary part, as required. 

In the incompressible case, equation (64) reduces to 


(67) 


dc 

dk 



(w - c) 2 


f dy 



and the imaginary part of dk/dc is 


( 68 ) 


Im( dk/dc) = Jr^ c |w c '^ 2 w c " ' J 0 2 dy 


(69) 
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NEUTRAL OSCILLATIONS 


From the discussions of the last section, another limitation is 
imposed on the occurrence of neutral subsonic disturbances . Since the 
wave speed of the neutral disturbance must he equal to the flow speed 
v s at the point where 


d / dw\ 
dy \ d yj 


= 0 


(70) 


subsonic disturbances cannot occur if the speed v s is supersonic rela- 
tive to either stream. Calculations of v s are made for the case of 
uni form total enthalpy. Instead of using equation (70) directly, it 
is found convenient to transform it into the form 


-1 

= 0 with ijr = f(£) and w = f'(£) (71) 


where \|r is the stream function, since the relation w(\|r) is identical 
with that in the incompressible case. The results of these calculations 
are shown in table IH. In comparing these values against condition (54), 
it is found that the condition 


& 

cfy‘ 


1 + 


(7 - liM-L 2 


- v2 


1 

/ 


M, 


< v. 


is never violated, hut the values of w s below the solid horizontal 
lines in the table are too high to satisfy the condition 


w c 



Thus, for a given Mach number in the stream M 2 , the Mach number 
in the other stream can become so high that the speed of flow corre- 
sponding to condition (70) becomes supersonic relative to the slower 
free stream. Subsonic disturbances then cannot exist. 
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The critical case of sonic disturbance is reached when 



Calculations for this case are tabulated in table IV. The correspo nding 
values of M]_ and Mg are plotted as curve B in figure 4 to mark the 
limit of stability. Only the region to the left of this curve can have 
subsonic disturbances. It is seen that the condition is more restrictive 
than that obtained from general considerations alone. In particular, in 
the case where one stream is at rest, the curve B shows that the flow 
becomes completely stable at a free-stream Mach number of 1-7? in con- 
trast to the value 2.5 given by curve A. Since they are, however, not 
very much different, it may be surmised that the exact distributions' of 
temperature and velocity in the mixing region may not be too important 
in determining the over -all stability characteristics. 

It may be recalled that the restriction to the case of constant 
total enthalpy can be immediately removed by considering a moving obser- 
ver. Thus, for all cases with Prandtl number equal to unity, it is nec- 
essary only to convert the values of the speeds U^, U 2 , and w s 

involved. Condition (70) is not modified by the reference to a moving 
observer. 

With the neutral wave speed thus determined, equation (55) can be 
integrated to give the amplitude of the oscillations. For this purpose, 
it is necessary to find the proper value for a. This can be done by 
several trials, with the first approximation given by the ratio of the 
integrals in equation (59) . The ratio will yield the characteristic 
value only when the function f (y) is the characteristic function, but 
it is known that any reasonable approximation to it will give a very 
close approximation to the characteristic value. 

Calculations of these neutral oscillations are carried out for 
several Mach numbers in one stream with the other at rest. These are 
given in table V and figure 6. 


SELF-EXCITED OSCILLATIONS 


The formulas given in the section "General Stability Characteristics 
in the Inviscid Case" have been used to calculate the characteristics of 
self-excited oscillations. It Is found that 


— = 0.177 - 0.2091 
da 
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for the case M^_ = 1 and 

— = 0.093 - 0.287i 
da 

for the incompressible case. The results are shown diagrammat i cally in 
figure 7. It is seen that the extent of amplification is fairly large. 
There is also indicated a decrease of amplification with increasing Mach 
number . 

It would be easy to calculate the amplitude function of these self- 
excited oscillations by using the characteristic values obtained above. 
This was not carried out because of limitations in time. 


NATURE OF OSCILLATIONS IN LIMIT OF 
DIFINITE REYNOLDS NUMBER 


Calculations in the section "Neutral Oscillations" indicate that 
the neutral oscillation has quite a simple amplitude distribution. In 
fact, it does not show any node. The amplified disturbances are expected 
to show similar characteristics. The mathematical theory (references 5 
and 6) , however, indicates that the damped oscillations behave in a much 
more complicated manner. It is concluded that there is always a finite 
viscous region in the interior of the fluid, no matter how large the 
Reynolds number may be. In fact, a minimum width of this region 1 b 
determined. This has to do with the crossing of the solid lines by the 
real axis of y in figure 5. The solution in the part of the real axis 
between the solid lines shows exponential behavior - and therefore viscous 
nature - if the solution in the outside parts shows the inviscid behavior. 
This type of conclusion has been reached in reference 1. However, it 
was possible only to suggest that such viscous behavior would occur at 
the solid lines. The improved theory shows that it must occur through- 
out the region in between. It is important to note that the complex con- 
jugate of an amplified solution does not represent a damped osc illa tion 
and vice versa, although this conclusion can be easily reached by a 
cursory exa min ation of the inviscid equations. The damped oscillations 
do not satisfy the inviscid equation al 1 along the real axis; otherwise, 
they could not take on the proper branch of the logarithm as specified 
in the section "General Study of Small Disturbances in a Nearly Parallel 
Flow Field in a Compressible Fluid." They exhibit a behavior very much 
like that of the vorticity distribution in fully developed turbulence 
flow. For large Reynolds numibers, there is one part of the space where 
the vorticity is highly concentrated; in another part, there is very 
little vorticity. This illustrates the two kinds of limiting behavior 
of a viscous fluid in the limit of infinite Reynolds n umb er; In one 
part of the field, the inviscid behavior is approached; in another part, 
it becomes highly oscillatory spatially. 
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CONCLUSIONS 


From a study of the stability of the mixing of two parallel streams 
in a gas, the following conclusions may he drawn: 

1. If the relative speed of the two parallel streams exceeds the sum 
of their velocities of sound, subsonic oscillations cannot occur, and 
the mixing region may be expected to be stable with respect to small 
disturbances . 

2. A further necessary condition for the possible occurrence of 
small subsonic disturbances is that somewhere in the field 



where y is a positional coordinate across the stream, p is the den- 
sity of the gas, and w is a dimensionless velocity distribution. This 
condition is usually satisfied for the present class of problems. 

3. If the speed of the flow at such a point is denoted by w B , 
then the field of flow can execute a neutral wavy oscillation having 
a finite wave length and propagating with the speed c = w s if and 
only if w s is subsonic relative to both streams. There is no other 

possible neutral oscillation. This leads to a more strict condition 
of stability than that given by conclusion 1 . 

I 4 .. Under the above conditions, the field of flow can execute 
amplified wavy oscillations having wave lengths longer than that of 
the neutral oscillation. Oscillations having shorter wave lengths axe 
damped. (This specific form of the conclusion was not obtained for 
the boundary layer at a solid surface in NACA TN 1115, but its validity 
can be shown by the present method.) The extent of amplification in 
such cases is fairly large. 

5. At large Reynolds numbers, the amplified disturbances are 
essentially free from the effect of viscosity. On the other hand, 
disturbances with finite dampi n g are expected to exhibit a highly 
oscillatory behavior over a finite region in the field of flow. This 
is similar to the structure of the vorticity field in fully developed 
turbulence . 

6 . For the case of constant enthalpy with one stream at rest, the 
wave length of the neutral disturbances increases with increasing Mach 
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number in the other stream. The flow becomes completely stable at a 
free -stream Mach number of 1.7. This is more restrictive than the 
value 2.5 obtained by applying conclusion 1 to the present case. 


Massachusetts Ins titute of Technology 
Cambridge, Mass., July 30, 1952 
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APPENDIX 

SYMBOLS 


The quantities hearing subscript 1 in the last column are the 
dimensional quantities in the first free stream. Corresponding quan- 
tities in the second stream hear a subscript 2. The quantities with- 
out a prime satisfy the equations of steady motion; those with primes 
satisfy the disturbance equations. 

Dimensional Dimensionless Reference 

quantities quantities quantities 

Positional coordinates: 

x* x l = 2\ 2V 1 x*/U 1 


y* 

y = V2tf2~ 

l 

Time: 

t* 

t 

Z/Ui 

Velocity components 

in directions of x- and y-axes. 

respectively: 

u* + u*' 

w(y)+f(y)e ia ( x - ct ) 

Ul 

V* y* * 

-a0(y) eLx(x-ct) 

Ul 

Density of gas: 

p* + p*» 

p(y)«(y>e ta(x - ot > 

Pi 

Pressure of gas: 



p* + p*» 

p(y)+«(y)e ia ( x - ct ) 

Pi 

Temperature of gas: 

T* 

T(y) +9 (y) e ia ( x " ct ) 

T 1 


Coefficient of viscosity of gas: 
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Dimensional Dimensionless 

quantities quantities 

Thermal conductivity: 

k* b/Pr 


Wave number of disturbance : 

a* = a = 2 it/\ 

Phase velocity of disturbance: 

c* c 


Specific heat at constant volume: 


Specific heat at constant pressure:- 
°P 7 


Gas constant per gram: 

R* 7-1 

Reynolds number 

R = PpUj_iy/ 

Mach number 

M 1 = Ui^7R*Ti 


Reference 

quantities 


C P P 1 


l - 1 


U 1 


°v 

°v 




Prandtl number 


Pr = CpP*/k* 
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TABLE I 

VELOCITY DISTRIBUTIONS IN INCOMPRESSIBLE CASE - Continued 


(d) For X = 0.8 and 0 = 0.860 


n 

V 

n 

V 


w 

-9.0 

0.111 

-2.6 

0.209 

1.4 

O.813 

-8.8 

.111 

-2.4 

.224 

1.6 

.842 

-8.4 

• 111 

-2.2 

.240 

1.8 

.869 

-8.0 

• 111 

-2.0 

.258 

2.0 

.893 

-7.6 

.112 

-1.8 

.279 

2.2 

.914 

-7.2 

.112 

-1.6 

.301 

2.4 

.932 

-6.8 

•113 

-1.4 

.326 

2.6 

• 947 

-6.4 

.115 

-1.2 

.352 

2.8 

.960 

-6.0 

■ 117 

-1.0 

.382 

3-0 

• 970 

-5.6 

.120 

-.8 

.413 

3-2 

• 977 

-5.2 

.124 

-.6 

.446 

3 -b 

.984 

-4.8 

.129 

-.4 

.481 

3.6 

.988 

-4.4 

.137 

-.2 

.518 

3-8 

• 992 

-4.0 

.146 

0 

.556 

4.0 

• 994 

-3-8 

.152 

.2 

.594 

4.2 

.996 

-3.6 

.159 

.4 

•633 

4.4 

.998 

- 3 -b 

.167 

.6 

.672 

4.6 

• 999 

-3-2 

.175 

.8 

.709 

4.8 

-999 

-3.0 

.185 

1.0 

.746 

5.0 

1.000 

-2.8 

.196 

1.2 

.780 
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TABLE II 

VELOCITY AND TEMPERATURE DISTRIBUTIONS 
IN COMPRESSIBLE CASE 


[One stream at rest; X = 1; 0 



(a) M]_ = 0.5 



V 

T 

n 

V 

T 

-12.6 

-12.2 

-u.8 

0 

1.050 

- 1.4 

-1.2 

-1.0 

0.261 

.289 

.320 

H H H 

b b b 

- 11.4 

-n.o 

-10.6 



-.8 

-.6 

-.4 

.352 

.386 

.423 

1.044 

1.043 

1.041 

-10.2 



-.2 

.461 

1.039 

-9.8 



0 

.500 

1.038 

- 9.4 

0 

1.050 

.2 

.540 

1.035 

-9.0 

.001 


.4 

.581 

1.033 

-8.6 

-8 2 

.001 

.002 

.003 


.6 

.8 

1.0 

.622 

.662 

.702 

1.031 

1.028 

1.025 

-7.8 


-7 .4 

.004 


1.2 

.740 

1.023 

- 7.0 

.006 


1.4 

.776 

1.020 

-6.6 

.009 


1.6 

.809 

1.017 

-6.2 

.013 


1.8 

.840 

1.015 

- 5.8 

.019 


2.0 

.868 

1.012 

- 5.4 

.025 


2.2 

.892 

1.010 

- 5.0 

• 033 

1.049 

2.4 

.914 

1.008 

- 4.6 

.043 


2.6 

.932 

1.007 

- 4.2 

.054 


2.8 

• 947 

1.005 

- 3-8 

.069 


3-0 

.960 

1.004 

- 3-4 

.086 


3-4 

• 978 

1.002 

- 3-0 

.108 


3-8 

.989 

1.001 

-2.8 

.121 


4.2 

.995 

1.001 

-2.6 

.136 


4.6 

• 998 

1.000 

- 2.4 

.152 


5.0 

• 999 

1.000 

-2.2 

.170 


5.4 

1.000 

1.000 

-2.0 

.190 

1.048 

5.6 

1.000 

1.000 

Ea 

.211 

.235 

1.048 

1.047 

6.0 

1.000 

1.000 



V 


T 


v 


T 


V 

T 


V 

0 

1.200 


0.282 

.309 



1 

.336 

.366 

.397 

.430 

0 

1.200 

-.2 

.465 

.001 


0 

.500 

.002 


.2 

.537 

.002 


.4 

.574 

.004 


.6 

.611 

.005 


.8 

.648 

.007 


1.0 

.685 

.010 


1.2 

.721 

.014 


1.4 

.756 

.019 


1.6 

.789 

.024 


1.8 

.820 

.031 


2.0 

.848 

.039 


2.2 

t- 

00 

• 

.048 


2.4 

.897 

.059 


2.6 

.918 

.072 

1.199 

2.8 

.935 

.088 

1.198 

3.0 

.950 

.108 

1.197 

■ 3.k 

• 971 

.132 

1.197 

3.8 

.985 

.145 

1.196 

4.2 

• 993 

.160 

1.195 

4.6 

• 997 

.177 

1.194 


.999 

.195 

1.192 


1.000 

.214 

1.191 



.235 

.258 

1.189 

1.187 

6.0 



1.184 

1.181 

1.177 

1.173 

1.168 ■ 

1.163 

1.157 

1.150 
1.142 ' 
1.134 
1.125 
1.116 
1.106 
1.096 
1.086 
1.076 
1.066 
1.056 
1.047 
1.039 
1.032 
1.025 
1.020 
1.011 
1.006 
1.003 
1.001 
1.001 
1.000 
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TABLE II 

VELOCITY AND TEMPERATURE DISTRIBUTIONS 


IN COMPRESSIBLE CASE - Concluded 


(c) M x = 1.5 


HH 

V 

T 

T) 

V 

T 

-12.6 

0 

1.450 

-1.4 

0.310 

1.407 

-12.2 

.001 


-1.2 

•334 

1.400 

-n.8 

.001 


-1.0 

.359 

1.392 

-n.4 

.002 


-.8 

.385 

1.384. 

-n.o 

.002 


-.6 

.412 

1.374 

-10.6 

.003 


-.4 

.440 

1.363 

-10.2 

.004 


-.2 

.470 

1.351 

-9.8 

.006 


0 

.500 

1.338. 

-9.4 

.007 


.2 

.531 

1.323 

-9.0 

.009 


.4 

.564 

1.307 

-8.6 

.012 


.6 

.596 

1.290 

-8.2 

.016 

— 

.8 

.629 

1.272 

-7.8 

.020 


1.0 

.662 

1.253 

-7-4 

.025 • 


1.2 

.694 

1.233 

-7.0 

.030 


1.4 

.727 

1.213 

-6.6 

• 037 

1.449 

1.6 

.758 

1.192 

-6.2 

.044 

1.449 

. 1.8 

.788 

1.171 

-5.8 

.052 

1.449 ' 

2.0 

£4 

CO 

• 

1.150 

-5.4 

.062 

1.448 

2.2 

.844 

1.130 

-5.0 

.073 

1.448 

2.4 

.869 

1.111 

-4.6 

.086 

1.447 

2.6 

.891 

1.093 

-4.2 

.102 

1.445 

2.8 

.911 

1.076 

-3.8 

.120 

1.444 

. 3-0 

.929 

1.062 

-3-4 

.142 

1.441 

3-4 

.957 

I.O38 

-3-0 

.167 

1A38 

3.8 

.976 

1.021 

-2.8 

.180 

1.435 

4.2 

.988 

1.011 

-2.6 

.196. 

1.433 

4.6 

.994 

1.005 

-2.4 

.212 

1.430 

5.0 

• 998 

1.002 

-2.2 

.229 

1.426 

5.4 

• 999 

1.001 

-2.0 

.248 

1.422 

5.6 

1.000 

1.000 

-1.8 

.267 

1.418 

6.0 

1.000 

1.000 

-1.6 

.288 

1.413 
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TABLE HI 

WAVE SPEEDS OF POSSIBLE NEUTRAL SUBSONIC DISTURBANCES 



v 8 


u*/Uj_ defined "by equation (70)J 



X 

0.2 

0.4 

0.6 

0.8 

1.0 

0 

0.8314- 

0.732 

0.657 

0.613 

0.576 

.5 

.845 

•743 

.678 

.632 

.594 

1 

.851 

.762 

.704 

.664 

.633 

2 

.878 

.823 

•796 

.772 

.762 

5 

.949 

.940 

.941 

.938 

.938 

10 

.984 

.981 

.981 

.982 

.982 


TABLE IV 

CONDITIONS ASSOCIATED WITH NEUTRAL SONIC DISTURBANCE 


Conditions 

X 

0.2 

0.4- 

0.6 

0.8 

1.0 

M 1 

>10 

4.0 

2.7 

2.0 

1-7 

“2 

>1.92 

.90 

.44 

.167 

0 

c 

«1 

.91 

.82 

.78 

.73 
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TABLE V 

WAVE NUMBER AND AMPLITUDE OF NEUTRAL OSCILLATIONS 

[y = n/2if2, e = i.i6cT] 

(a) Wave number of neutral oscillations 
for various Mach numbers 


% 

Wave number, cud 

0 

0.k59 

.5 

.k5k 

1.0 

•37^ 

1.5 

.324 
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TABLE V 

WAVE NUMBER AND AMPLITUDE OF NEUTRAL OSCILLATIONS - Continued 


(t>) Amplitude of neutral osc illa tions 
for M-j- = 0 
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TABLE V 

WAVE NUMBER AND AMPLITUDE OF NEUTRAL OSCILLATIONS - Concluded 

(c) Amplitude of neutral oscillations 
for = 1 


y 

0 

d0/dy 

y 

0 

d0/dy 

y 

0 

d0/dy 

- l.k 

0.003 

0.003 

-0.3 

0.551 

0.284 

2.9 

0.122 

-0.106 

-7-0 

.005 

.003 

-.2 

• 579 

.273 

3.0 

.112 

-.098 

-6.6 

.006 

.005 

-.1 

.605 

.254 

3-i 

.103 

-.090 

-6.2 

.008 

.006 

0 

.629 

.227 

3.2 

.094 

-.082 

-5.8 

.011 

.008 

.1 

.650 

.190 

3-3 

.087 

-.076 

-5.lt- 

.015 

.011 

.2 

.667 

.142 

3-4 

.079 

. -.070 

-5.0 

.020 

.015 

.3 

.678 

.085 

3.6 

.067 

-.059 

-4.6 

.026 

.020 

.4 

.684 

.020 

3.8 

.056 

-.049 

-4.2 

.035 

.026 

.5 

.682 

-.050 

4.0 

.048 

-.042 

-3-8 

.048 

.035 

.6 

.674 

-.121 

4.2 

.040 

. -.035 

-3-4 

.064 

.047 

.7 

.658 

-.189 

4.4 

•.034 

-.030 

-3-0 

.086 

.063 

.8 

.636 

-.249 

4.6 

.029 

-.025 

-2.6 

.115 

.084 

• 9 

.609 

-.297 

4.8 

.024 

-.021 

-2.2 

.153 

.112 

1.0 

.577 

-•331 

5.0 

.020 

-.018 

-2.0 

.177° 

.129 

1.1 

.543 

--351 

5.2 

.017 

-.015 

-1.9 

.191 

.138 

1.2 

.507 

-.359 

5.4 

.014 

-.013 

-1.8 

.205 

.148 

1-3 

.471 

-.35 6 

5.6 

.012 

-.011 

-1.7 

.220 

.158 

1.4 

.436 

-.344 

5.8 

.010 

-.009 

-1.6 

.237 

.169 

1.5 

.403 

-.328 

6.0 

.009 

-.008 

-1.5 

.254 

.180 

1.6 

• 371 

-.308 

6.2 

.007 

-.006 

-1.4 

.273 

.192 

1.7 

.341 

-.288 




-1.3 

.292 

.204 

1.8 

•3-3 

-.267 




-1.2 

•313 

.217 

1.9 

.288 

-.246 




-1.1 

.338 

.229 

2.0 

.264 

-.227 




-1.0 

.359 

.242 

2.1 

.242 

-.209 




-.9 

.384 

.254 

2.2 

.222 

-.192 




-.8 

.410 

.265 

2.3 

.204 

-.177 




-.7 

• 437 

.275 

2.4 

.187 

-.163 




-.6 

.465 

.282 

2.5 

.171 

-.150 




-.5 

.493 

.287 

2.6 

.157 

-.137 




-.4 

.522 

.288 

2.7 

.144 

-.126 







2.8 

.132 

-.116 




















(a) With "both streams moving. (h) Wd,th one stream at rest. 

Mi and M2; Mach numbers Ml, Mach number of moving 

of two streams. stream; M2 = 0; 7 = 1.40. 


Figure 3.- Stable zones for various temperatures and Mach numbers, 
Tg/Ti, ratio of absolute temperatures of two streams. 
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(c) For Mi = 1.0 


and a0 = O. 37 A. 



(d) For Mi =■ 1.5 and a 9 = 0.32k. 
Figure 6.- Concluded. 
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(b) For Mi = 1. 

Figure 7 *- Extent of amplification at infinite Reynolds n umb er , 

a, .wave number; — , wave length in y-units; y = -i- r) = — ; 

2 1/2 2l/2 /v^/Uy 

Rj Reynolds number; e, dimensionless complex wave speed ^c^ + ic^). 
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